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We have observed Bragg scattering of photons from quantum degenerate 87 Rb atoms in a three- 
dimensional optical lattice. Bragg scattered light directly probes the microscopic crystal structure 
and atomic wavefunction whose position and momentum width is Heisenberg-limited. The spatial 
coherence of the wavefunction leads to revivals in the Bragg scattered light due to the atomic Talbot 
effect. The decay of revivals across the superfluid to Mott insulator transition indicates the loss of 
superfluid coherence. 
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Ultracold atomic gases are an ideal system in which 
to study many-body phenomena because of the relative 
ease with which parameters in the model Hamiltonian 
can be tuned across a wide range [U E]. Such studies 
have resulted in a better understanding of various phase 
transitions such as the Berezinskii-Kosterlitz-Thouless 
transition in two dimensional systems [3], the BEC-BCS 
crossover of interacting fermions [4] and the superfluid 
to Mott insulator transition in a three-dimensional lat- 
tice [5]. One major goal of this field is to realize spin 
phases such as antiferromagnetic states to explore quan- 
tum magnetism and its interplay with high-temperature 
superconductivity [Sj. 

Concurrently, there are numerous efforts to develop 
techniques to probe and understand the atomic ensem- 
ble once a new type of ordering is achieved. One recent 
development is the realization of single-site resolution of 
atoms in two-dimensional optical lattices [3 |S] • An al- 
ternative method to measure in situ spatial ordering is 
the technique of Bragg scattering, often used in a con- 
densed matter context to determine crystal structure [9] . 
In particular, Bragg scattering with neutrons led to the 
discovery of antiferromagnetism in the solid state [TU] 
and with x-rays led to the discovery of the double helix 
structure of DNA [TT] . 

Bragg scattering relies on the interference of waves 
scattered coherently from an ensemble of particles. In 
particular when atoms are arranged in a periodic pat- 
tern in three spatial dimensions, there are certain angles 
of the incoming and reflected beams where scattering is 
dramatically enhanced compared to other angles. This 
has allowed crystallographers to use x-rays to determine 
the properties of crystals such as lattice geometry at the 
angstrom scale. We have applied this technique to ul- 
tracold atoms in a three-dimensional optical lattice by 
scattering near-infrared light where the atoms are spaced 
almost 10 4 times farther apart than those in typical con- 
densed matter samples. 

Pioneering works on Bragg scattering from cold atoms 
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FIG. 1. Bragg scattering as a probe of the spatial wavefunc- 
tion width. Bragg scattered intensity vs. lattice depth in 
units of the recoil energy. The right axis gives the correspond- 
ing root-mean-square width of the wavefunction squared. The 
solid line is a no-free-parameter curve given by the Debye- 
Waller factor. Error bars are statistical. Inset is a schematic 
of the setup for Bragg scattering. 



in optical lattices were done by the Hansch and Phillips 
groups using laser cooled atoms [12]. These lattices 
were sparsely populated and the atoms occupied several 
bands. In this Letter, we have used Bragg scattering 
to study bosonic atoms cooled to quantum degeneracy 
and placed in a three dimensional cubic lattice where the 
atoms occupy the lowest band. In particular, we have 
measured directly the Heisenberg-limited width of both 
position and momentum of the single ground state atomic 
wavefunction in an optical lattice. Furthermore, there is 
a revival of Bragg scattered light some time after releas- 
ing the atoms from the optical lattice, analogous to the 
optical Talbot effect. This signal gives a direct measure 
of the coherence of the superfluid component in the lat- 
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tice. 

The experimental apparatus has been described in de- 
tail elsewhere [13]. Briefly, laser cooling and evaporative 
cooling are used to achieve quantum degeneracy of 87 Rb 
atoms in the \F = 2,m F = —2) state which are loaded 
into a crossed dipole trap whose trap frequencies range 
between 30 and 160 Hz. Once quantum degeneracy is 
achieved, the optical lattices generated by a single 1064 
nm fiber laser are adiabatically ramped on. The lattices 
are calibrated by applying a 12.5 /is pulse and measuring 
the Kapitza-Dirac diffraction of the atoms and compar- 
ing this to theory. The system typically contains about 
10 5 atoms, which leads to up to 5 atoms per lattice site. 
The Bragg reflected light is detected on a CCD camera 
which images along a direction which satisifies the Bragg 
condition. 

The Bragg scattering condition for a three-dimensional 
cubic lattice is given by 2d sin 8 = n\ p where d is the 
spacing between lattice planes, is the angle between the 
incoming beam and the lattice planes, n is any positive 
integer and X p is the wavelength of the probe beam. For 
our experiment X p is 780 nm and d is 532 nm. With these 
conditions the only allowed angle 9 is 47° where n is one 
and corresponds to Bragg reflection off the [1 0] plane 
or any equivalent plane. A schematic of the probe beam 
with respect to the atoms is shown in the inset of Fig. [T] 
Since the full angular width of the Bragg scattered light 
is small (measured to be 4.1 ± 0.4°), a precise alignment 
of the incident beam had to be performed at the 1/min 
repetition rate of the experiment. In contrast, in one 
and two dimensions, diffractive light scattering occurs at 
any angle of incidence, as recently shown with atoms in 
a two-dimensional optical lattice [T3] . 

The probe beam at the Bragg angle had a typical power 
of 0.3 mW and beam diameter of 300 /im, large enough 
to illuminate the entire atomic cloud, and was detuned 
from the 5 2 S 1/2 ,F = 2 5 2 P 3/2 ,F' = 3 cycling tran- 
sition of 87 Rb by a few natural linewidths, where the 
natural linewidth T is 6 MHz. The detuning needs to 
be sufficient so that the light traverses the entire atom 
cloud. Although the absolute signal varies with the de- 
tuning, the conclusions reached in this Letter appear not 
to depend strongly on the amplitude or sign of the detun- 
ing. The probe beam was applied for a few microseconds, 
enough to obtain a signal but short enough to avoid heat- 
ing effects. 

As the lattice is increased the Bragg reflected signal in- 
creases as expected from the crystal ordering and tighter 
localization of the atoms (Fig. [lj. For this data, the 
lattice in the z direction, which is also the direction of 
gravity, was ramped to 15.Br, where Er — h 2 /(2mA 2 ) 
is the recoil energy, h is the Planck constant, m is the 
mass of 87 Rb and Al = 1064 nm is the lattice laser wave- 
length. Simultaneously the lattices in the horizontal di- 
rections were ramped to various lattice depths. Note that 
for our geometry the Bragg reflected intensity is insensi- 
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FIG. 2. Bragg scattering as a probe of the momentum 
wavefunction width. Bragg scattered intensity vs. the free- 
expansion time of the atoms after rapidly turning off the lat- 
tices from 15i?R. The decay in signal indicates a melting 
of the crystal structure, or in other words spreading of the 
atomic wavepacket with a given momentum uncertainty. The 
gray area is a no-free-parameter curve using the Debye- Waller 
factor taking into account the probe pulse duration of 5 fis. 
Error bars are statistical. 



tive to the lattice depth in the z direction because Bragg 
scattering occurs only in the horizontal xy plane, as we 
discuss later. 

The scattering of light by a collection of atoms can 
be modeled in the following way. In the limit of low 
probe intensity and low optical density of the cloud, 
the Born approximation can be used. Then the scat- 
tering cross-section da/d£l can be written as a product 
of one-dimensional Debye- Waller factors where da/dCl oc 
Y\ i exp(—(Ar i ) 2 Kf/2), where i = x,y, z is the index for 
the three dimensions, Ki is the magnitude of the i-th 
component of K where K = kj n — k out is the difference 
between the incoming probe beam and Bragg scattered 
wavevectors, and Ar^ = y/h/(mu)i) is the harmonic oscil- 
lator width where Ti is the reduced Planck constant and Ui 
is the trap frequency in the z-th direction, which depends 
on the optical lattice depth. The atoms are approximated 
as gaussian wavefunctions, which is a good approxima- 
tion for sufficiently large lattice depths. However, we 
find that even for low lattice depths where significant 
superfluid components are expected, this approximation 
describes the Bragg scattered signal well. 

For our experimental conditions, K z — because 
Bragg scattering occurs in the horizontal xy plane. The 
lattice depths are controlled in a way such that they are 
the same in both the x and y directions which leads to 
Ax = Ar x = Ar y . This simplifies the scattering cross- 
section to da/dfl ex exp(—(Ax) 2 K 2 /2). Thus Bragg 
scattering allows the study of the spatial extent of the 
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atomic wavef unction. 

The data can be compared to a no-free-parameter the- 
oretical line, assuming non-interacting atoms. The Bragg 
scattered intensity B(Nl) as a function of lattice depth 
N L is proportional to exp(- \ 2 L K 2 / (8tt 2 y^)) . Both 
the harmonic oscillator width Ax, which depends on 
the mass of the rubidium atom and trap frequency, and 
change in wavevector K are known. This theoretical line 
and the data is shown in Fig. [T] and we see good agree- 
ment. Thus the Bragg scattered light as a function of 
lattice depth can be well-described by a model where the 
atoms are assumed to be non-interacting gaussian atomic 
wavefunctions whose spatial width is determined by the 
lattice depth. 

Bragg scattering can also be used to probe the momen- 
tum width of the atomic wavefunctions. This is done by 
measuring Bragg reflection as a function of the expan- 
sion time of the atoms between a rapid lattice turn off 
and Bragg probe. In particular, the two horizontal lat- 
tices were turned off from 15.Br to 0£r in less than 1 
fxs and the lattice in the z direction was kept at 15Ar. 
Turning off the lattices allows the atomic wavepackets in 
each individual lattice well to expand freely in those di- 
rections. The data in Fig. [2] shows that the signal has 
decayed in 30 /is. The time it takes to lose crystal struc- 
ture is roughly the time it takes for the atoms to move 
over half of a lattice distance. 

The Debye- Waller factor can be used to determine 
more precisely how the Bragg reflection should behave 
as a function of the expansion time with no free pa- 
rameters, assuming gaussian atomic wavepackets. This 
makes use of the well-known time-dependent behavior of 
a Heisenberg-limited gaussian wavepacket in free space 
which can be written as (Ax(t)) 2 = (Ax) 2 + (Ap) 2 t 2 /m 2 
where Ax and Ap are the uncertainty of position and 
momentum respectively at the initial time and t is the 
expansion time [15 . In a previous paper where the atoms 
were not quantum degenerate, the momentum distribu- 
tion was assumed to be determined by temperature [TJ] . 
The results in this Letter show that the momentum un- 
certainty is Heisenberg-limited, meaning (Ax) 2 (Ap) 2 — 
h 2 /4, so the Bragg scattered light B(t) as a function of 
expansion time decays as exp(— (Ap) 2 K 2 t 2 / (2m 2 )). This 
curve is also shown in Fig. [2] The curve is broadened by 
taking into account the probe beam duration which was 
5 fxs and shows good agreement with data. This anal- 
ysis shows that by releasing the atoms from a lattice, 
we can directly probe the in situ time evolution of the 
ground state atomic wavefunctions with spatial extent of 
tens of nanometers. Furthermore, the atomic spatial and 
momentum widths are seen to be Heisenberg-limited. 

Coherent many-body non-equilibrium dynamics can 
also be studied with Bragg scattering after a sudden re- 
lease of the atoms from the optical lattice in three dimen- 
sions. At low lattice depths the Bragg scattered signal 
shows revivals as a function of expansion time because of 
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FIG. 3. Revivals of Bragg scattered light. Bragg scattered 
intensity vs. the expansion time in three dimensions for three 
different initial lattice depths: 5-Er (circles), 8En (triangles), 
and 15-Br (squares). Each data set is normalized to the Bragg 
intensity at the initial time. Revivals at lower lattice depths 
indicate coherence of the atoms across multiple lattice sites. 
The lines are phenomenological fits to exponentially decaying 
sine waves. Different lattice depth data are offset for clarity 
and the error bar is represenative. Inset is a solution to the 
one-dimensional Gross-Pitaevskii equation with experimental 
parameters and an initial state of a chain of gaussian wave- 
functions, showing the revivals of the density distribution as a 
function of time. The white dotted line represents the revival 
time. 



the rephasing of the superfluid atomic wavefunctions that 
were originally confined in the lattices. Therefore, these 
revival signals should give a measure of the superfluid 
order parameter. Furthermore, the revivals are analo- 
gous to the optical Talbot effect, whose atomic version 
was observed previously for a thermal beam |16j and for 
a quantum degenerate gas in the temporal domain |17j . 
Here we observe the atomic Talbot effect for an interact- 
ing quantum degenerate gas in a three-dimensional op- 
tical lattice. Revivals can be seen in Fig. [3] and become 
less pronounced as the lattice depth increases. 

In particular, the characteristic revival time is deter- 
mined by t = h j (ft 2 (2k) 2 j (2m)) and is 123 fj,s for our 
parameters, where 2k corresponds to the wavevector of 
the matter wave. However, with Bragg scattering we ex- 
pect the first revival at half that time of 61 /its, which is 
what we observe in Fig. [3] This can be understood by 
realizing that the atomic wavepackets need only travel 
half the lattice spacing to constructively interfere with 
the wavepackets traveling in the opposite directions from 
the nearest neighbor sites. 

The superfluid coherence as a function of lattice depth 
can be studied by comparing the Bragg reflected signal 
when the atoms are in the optical lattice to the signal 
at the first revival. The revival signal as a function of 
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FIG. 4. Bragg scattering revival as a probe of superfluid co- 
herence. Top: Bragg scattered revival signal at 60 /is vs. 
lattice depth. Bottom: The top plot normalized to the Bragg 
scattered signal without expansion. We see a decrease in re- 
vival, and consequently superfluid coherence, as the lattice is 
increased. Error bars are statistical. 



lattice depth is shown in Fig. |4j where the superfluid 
to Mott insulator transition is expected to occur around 
13-Er, [5]. In a non- interacting system without any dissi- 
pation, one would expect complete revivals. The decrease 
in revival fraction as the lattice depth increases is consis- 
tent with the picture that as interactions increase, phase 
fluctuations among neighboring lattice sites increase and 
consequently the superfluid fraction decreases. In the 
Mott insulating phase, revivals should completely disap- 
pear, except for a very weak signal at intermediate lattice 
depths due to particle-hole correlations |18j . 

To understand the measured revival decay, we have nu- 
merically studied the one-dimensional Gross-Pitaevskii 
equation that assumes interacting matter waves at zero 
temperature. The simulations show that interactions 
and finite size effects have negligible effect on the de- 
cay of revivals. Empirically, randomizing the phase be- 
tween neighboring sites reduces the revival fraction. The 
loss of phase coherence across lattice sites could be due 
to factors such as temperature, beyond Gross-Pitaevskii 
equation effects, or technical factors. We have considered 
quantum depletion [TS] , but the calculated depletion frac- 
tion at 5-Er is 2%, too low to account for the observed 
decay. Future work should provide a more complete pic- 
ture of the decay of revivals. 

Note that Bragg scattered revivals are complementary 
to the observation of diffraction peaks in time-of-flight 



absorption images [5]. Both are based on matter- wave 
interference due to superfluid coherence: revivals at short 
expansion times and diffraction at long expansion times. 
Diffraction is a far-field effect only possible for finite 
size samples, whereas revivals are a near-field bulk ef- 
fect. The order of the revival or the angular resolution of 
the diffraction peaks determine whether these techniques 
probe short-range or long-range spatial coherence. Fur- 
ther studies are needed to determine which technique is 
more sensitive to certain aspects of the superfluid to Mott 
insulator transition. 

In this Letter we have not focused on the effects of oc- 
cupation number on Bragg scattering. In principle the 
atomic wavefunctions are more extended for higher oc- 
cupation numbers, but this effect is small for our pa- 
rameters. However, light scattering at higher occupation 
numbers will have an inelastic component due to collid- 
ing atoms or photoassociation (2U]. The dependence of 
Bragg scattering on lattice depth suggests that these ef- 
fects are not dominant. 

After Bragg scattering has been established as a probe 
for Mott insulators, it can now be applied to study other 
types of quantum phases such as antiferromagnetic or- 
dering in both the occupation number and spin sec- 
tors [2lJ [22]. Although temperatures required to real- 
ize spin ordering are on the order of tens to hundreds 
of picokelvins, recent experimental results suggest a way 
forward to lower the temperature of a two-component 
Mott insulator |25j . 

In conclusion, we have observed Bragg scattering of 
near-resonant photons from a quantum degenerate Bose 
gas in a three-dimensional optical lattice. We have shown 
that this technique probes not only the periodic structure 
of the atoms, but also the spatial and momentum width 
of the localized atomic wavefunctions and the superfluid 
coherence. 
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